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Abstract 

This short note is based on an early draft of a paper we wrote describing our 
measurement of the hyperfine quenching rate of the clock transition in Yb-171, see 
Phys. Rev. Lett. 113 033003 (2014) Hj. We discuss the El and Ml spontaneous decay 
rates of the an emitter residing inside of a real cavity carved out of a vast, uniform, 
homogenous, isotropic, linear, lossless, dispersionless, and continuous medium. The 
ratio of the medium rate to vacuum rate is given by Tm/To = [G{u)]‘^n^/u, where 
G{u) = 3m/(2m + 1) is the local field correction factor, n = ye/i/(eo/io) is the index of 
refraction of the medium, e(eo) is the electric permitivity of the medium (vacuum), 
p(/io) is the magnetic permeability of the medium (vacuum), and m = e/eo for El 
transitions or m = po/ii for Ml transitions. 
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1 Introduction 

Spontaneous emission occurs when an emitter, such as an excited atom or molecule, cou¬ 
ples to the zero-point energy fluctuations of the electromagnetic field and undergoes a 
quantum transition resulting in a photon. The rate of these transitions, or equivalently 
the lifetime of the excited state, depends not only on the intrinsic properties of the emitter, 
but also on the properties of the surrounding environment. Since Purcell first pointed this 
out in 1946 |[2|, there have been many experimental demonstrations of this so-called "Pur¬ 
cell effect" m m 13 and it is one of the hallmarks of cavity quantum electrodynamics. In 
recent years, materials are being engineered with a microscopic structure designed specif¬ 
ically to enhance the spontaneous emission rate in order to realize, among other things, 
miniature lasers and single photon sources HUZilSll- Even with all this prodigious effort, 
there is still considerable tension in the literature about how the spontaneous emission 
rate of an emitter is modified by the medium that surrounds it |[9|. 
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2 Previous Work on Dielectric Media 


This problem has been most extensively studied for electric dipole (El) transitions in 
homogenous, uniform, isotropic, lossless, and linear dielectric media with dielectric con¬ 
stant e and index of refraction n = \Je/eo. The main difficulties in addressing this question 
are understanding the way in which the emitter modifies the electromagnetic field in the 
surrounding medium and the way in which the medium modifies the properties of the 
emitter. The former being the main challenge in choosing the appropriate theoretical 
framework within which to calculate the spontaneous emission rate and the latter being 
the main challenge in interpreting experiments. 

The two main types of theoretical approaches are categorized as macroscopic and mi¬ 
croscopic. The main difference between these two approaches has to do with whether the 
medium is treated as continuous (macroscopic) or discrete (microscopic). In macroscopic 
theories, first, the electromagnetic field is quantized in the usual way llTOl ITTlIT^ [T3l , but 
now accounting for the proper normalization of the electromagnetic field fluctuations 
due to the permittivity of the medium. Then, to calculate the spontaneous decay rate, 
one applies Fermi's golden rule noting that the photon density of states scales as the 
wave number cubed k^, which implies an scaling. Putting this altogether, accounting 
for the medium modification of the emitter fluorescence spectrum UTIB . and introducing 
the local field correction IIT^ gives the ratio of the spontaneous decay rate in medium 
(qmX'm) to that in vacuum (To) 
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where u is the frequency of the transition, D is dipole matrix element, and the subscript 
m (0) refers to the medium (vacuum). The factor q^, known as the quantum efficiency of 
the transition, is the ratio of the spontaneous emission rate to the total excited state decay 
rate in medium. It is included to account for the possible presence of nonradiative decay 
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mechanisms induced by the medium. 

The factor G{u), known as the local field correction, is the ratio of the local field im¬ 
mediately surrounding the emitter and the macroscopic field in the dielectric medium far 
from the emitter. It is usually calculated by assuming the emitter is a point dipole that 
is located at the center of a dielectric sphere surrounded by the medium. Depending on 
dielectric constant of this sphere, G{u) is either Grcm = Su/{2u + 1) assuming the real cav¬ 
ity model (RCM) or Gvcm = (m -I- 2)/3 assuming the virtual cavity model (VCM). The choice 
of model is determined by how similar the emitter is to the medium, or, alternatively, by 
whether the emitter is located at a substitutional site (RCM) or an interstitial site (VCM) 
within a discrete lattice ||T6l . 

In microscopic theories, the medium atoms and emitter atoms are considered sep¬ 
arately and couple only via the electromagnetic field. These calculations typically rely 
on quantum many body calculations which attempt to account for correlations between 
nearby medium atoms induced by the presence of the emitter atom IfTTl IT^ [T^ l20l 1211 . 
The modification of the spontaneous decay rate in these theories is usually written in one 
of the two following ways 
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1 + CiaoN + C2{aQNY + 0{aoN)^ 

( 2 ) 


where A; = 0 or 1, now G is linear and not quadratic, oq is the electric polarizability of 
the medium atoms, N is the number density of the medium atoms, and ci & C 2 depend 
on the details of the calculation. We note that the units of uq are such that it satisfies the 
Clausius-Mossoti relation given by a^N = 3(n^ — l)/(n^ -|- 2). In all cases, a larger index 
of refraction is expected to result in a faster rate. 

The index of refraction dependence of the spontaneous decay rate has been measured 
in rare earth ion-complexes in liquids 1221 . dye molecules dissolved in water droplets 
suspended in hydrophobic liquids l23l UM . rare earth ion-complexes in dense gases 1251 , 
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nanoparticles in air & liquids 11261 , rare earth ions in binary glasses 1271128l 1291 , quan¬ 
tum dots in liquids l30l , and rare earth ions in vacuum and crystalline solids IMl |321 . 
Challenges in interpreting these results include determining the index of refraction of the 
composite systems, the quantum efficiency of fhe emiffer, and fhe changes in the dipole 
matrix element of the emitter due to the surrounding environment. In short, there is not 
a strong consensus in either the theory or experimental community about the correct so¬ 
lution to this problem. The current status of this problem has been recently reviewed by 
Dolgaleva & Boyd |33| . 

3 Simplifying Assumptions 

In order to make an unambiguous calculation, we restrict ourselves to El and Ml spon¬ 
taneous decay rates of fhe an emitter residing inside of a cavify carved ouf of a vast, uni¬ 
form, homogenous, isofropic, linear, lossless, dispersionless, and continuous medium. 
We assume that the emitter density is sufficiently dilute that we do not need to worry 
about the effect of a neighbhoring emiffer. A "vasf" medium simply means thaf the emit¬ 
ter is far from the boundaries of the medium. The constraint that the medium is uniform, 
homogenous, isofropic, linear, lossless, and dispersionless means fhaf fhe elecfromag- 
netic properties of fhe medium can be described by fwo scalar, real, frequency indepen- 
danf consfanfs: fhe elecfric permittibify e and fhe magnetic permeabilify p. The index 
of refraction of the medium is therefore n = y^ep/^eopo), where the subscript 0 labels the 
vacuum values. By placing the emitter inside of a cavity that is carved out of a contin¬ 
uous medium, we are able to unambiguously determine the local field correction facfor, 
which is given by fhe real cavify model as G{u) = 3u/{2u -1-1). In addition, we will simply 
assume that the oscillator strength and frequency of fhe fransitions are nof modified by 
the medium. 
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4 Atomic Electric Dipole Transitions 

We start by Fermi's Golden Rule which gives the tranistion rate due to an operator "H 
averaged over initial states and summed over final states: 

where h is the Planck constant divided by 27r, a{b) labels an initial (final) state, Na(b) is 
the total number of initial (final) states, and p is density of available states (number of 
states per unity energy). The transition matrix element for electromagnetic transitions in¬ 
volves the dot product between a dipole moment operator and an electromagnetic field. 
The dipole momentum operator acts only on the wave function representing the emit¬ 
ter, while the matrix element involving the electromagnetic field is related to the RMS 
fluctuations of the zero-photon states of the electromagnetic field. 

We'll first consider El transitons for atomic emitters where the electric dipole moment 
is d and the electric field is E: 
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where we've taken advantage of the orthogonality of electric dipole moment vector op¬ 
erator and the isotropy of space. 

In order to determine the matrix element related to the electric field, we first recall that 
the energy density of a classical electric field is given by: 


D-E 

2 



( 8 ) 


6 



where e is the dielectric constant. This energy density can be related to the zero point 
energy hum/‘2 of a simple harmonic oscillator summed over each mode labeled by m Il34l : 



where the 2 in the denominator indicates that the energy is partioned equally between the 
electric & magnetic fields, and V is the volume of the space occupied by all the modes. 
This sum and volume only make sense when combined with the density of photon states, 
which will be discussed shortly. What we've calculated so far involves the fluctuations 
of the electric field far from the emitter. However, as mentioned before, the electric field 
very close to the emitter are modified by a local field correction factor G = Eiocai/Ei^r. In 
our case, it is appropriate to invoke the real cavity model Il35l : 

G(«) = GrcmM = ^ (10) 

We note that G{u) varies from 1 to 3/2 for u>l. 

We now return to the density of photon states. First, the total number of photon states 
in the volume V with can contain a maximum wavelength A or alternatively minimum 
wave number k is given by: 

= (f)(1)' 

where the 2 indicates two polarization modes and we're assuming that V is very large. 
The relationship between wave number and frequency is given by: 

.A = (2..) (^) = j = / (12) 

where c is the speed of light in vacuum. The density of modes near k is given by Il36l : 
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Finally, we now stipulate that the only modes that contribute to the sum the ones 

with energy nearly equal \o \Ei, — Ea\ = hu = fikc/n. Putting this altogether gives: 
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where the final answer no longer depends on V. 

Before writing down the El spontaneous decay rate in medium, it is convenient to 
recall the following equations for the fine structure constant and the oscialltor strength: 
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where e is the elementary charge, m is the electron mass, and d = er. With these substitu¬ 
tions, we find the El spontaneous decay rate for an excited atom is: 
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5 Nuclear Spin Magnetic Dipole Transitions 

We can perform an analgous calculation for Ml transitions, which involves the dot prod¬ 
uct jl ■ B, where p is the magnetic dipole moment and B is the magnetic field. Following 
the energy density argument from before IIMIl , we find: 
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wher /i is the magnetic permeability of the medium. The local field correction factor has 
the same form with e —?• l//i and we find: 
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where po is the permeability of the vacuum. The product of the local 5-field matrix 
element squared and the photon density of state (which is the same) is 
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Putting this altogeher gives: 
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where G is now a function of po//^ arid the factor of 3 comes from the fact that the fluctu¬ 
ations of the magnetic field are isotropic. 

What is remaning now is the calculation of the magnetic dipole matrix element. For 
this calculation, we will assume that our emitter is a unpolarized nucleus with spin S 
and a magnetic moment given by |/2| = gpnS, where pat is the nuclear magneton. The Ml 
transitions that we considering are between the Zeeman levels of the spin that are split 
by some bias magnetic field giving energy levels hum = gpN'mB. WeTl start by writing 
the magnetic dipole vector operator in terms of spin ladder operators: 


jl = gpnS 
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Without loss of generality, we will assume the higher energy state is |m — 1), while the 
lower energy state is \m), which implies that (7 > 0. By noting that (m — 1| 5*+ \m) = 
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{m — l\Sz\m) = 0, the matrix element is then given by: 


(m — 1| Jl\m) 
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(24) 


By noting that (m — 1| S'. \m) = \JS{S + 1) — m{m — 1), = S'(S' + 1)(2S' + l)/6, 

and J2m=-s ^ = 0/ we find that the amplitude squared averaged over all initial values is: 
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Putting this altoghether gives the Ml spontaneous decay rate for an unpolarized nuclear 
spin S' is a magnetic field B\ 
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For a sense of scale, using the gf-factor of an electron in nuclear magneton units \ge\= 3681, 
i? = 10 T, and p = po/ we find r^\ = 1.43 yrs. If we were to put in ^ 1, then the resulting 
decay lifetime would be several orders of magnitude longer than the age of the Universe. 


6 Conclusion 

Determining the spontaneous decay rate in medium is complicated. The medium can: 

1. shift the emitter's transition frequency 

2. change the emitter's transition matrix element 

3. introduce non-radiative quenching mechanisms 
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The first two directly alter the spontaneous decay rate. The third makes the interpretation 
ot experimental data more challenging. Furthermore, the medium alters the effect of the 
fluctuations of the zero-photon state of the electromagnetic field by: 


1. increasing the density of photon states by a factor of 

2. altering the magnitude of the field fluctuations far from the emitter by e~^ or p 

3. amplifying the magnitude of the field flucuations near the emitter by a local field 
correction factor G{u) where u = e/or u = hq/^ 

The first one is well known. The last two is a matter of considerable debate in the litera¬ 
ture. The situation becomes even more complicated if the medium is: 

1. small compared to the transition wavelength: edge effects due boundary conditions 

2. nonuniform: e and p depend on position 

3. inhomogenous: e and p have to be averaged over different medium species 

4. anisotropic: e and p depend on the direction of emmission 

5. nonlinear 

6. lossy: e and p are complex which implies the medium can reabsorb the emitted light 

7. dispersive: e and p depend on the transition frequency 

8. discrete: some ambiguity about how to quantize the electromagnetic field 

If the emitter can be thought ot residing in a real cavity carved out ot a simple continuous 
medium and the emitter is unaffected by the medium, then the ratio of the in-medium 
decay rate to the vacuum decay rate is: 
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where n = ye/x/(eo/io) is the index of refraction of the medium and u = e/e^ for El fran- 
sitions or M = /io//i for Ml fransitions. For the special case when /i = hq, then n = y^e/eo, 
u = for El transitions, and u = 1 for Ml transitions, and medium to vacuum decay 
rate ratios are now: 
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which agrees with IITOl [22j but not IITlI . Based on the energy density scaling argument 
used in 113411 , we believe our solution is correct. 


References 

[1] C.-Y. Xu, J. Singh, J. C. Zappala, K. G. Bailey, M. R. Dietrich, J. P. 

Greene, W. Jiang, N. D. Lemke, Z.-T. Lu, P. Mueller, and T. P O'Goimor. 

Measurement of the Hyperfine Quenching Rate of the Glock Transition in Yb, 

Phys. Rev. Lett., 113:033003, Jul 2014. 

[2] E. M. Purcell. BIO. Spontaneous Emission Probabilities at Radio Frequencies. Phys. 

Rev., 69:681, June 1946. 

[3] K.H. Drexhage. Influence of a Dielecfric Interface on Fluorescence Decay Time. Jour¬ 
nal of Luminescence, 1-2:693-701,1970. 

[4] E. Yablonovifch, T. J. Gmitter, and R. Bhat. Inhibited and Enhanced Spontaneous Emission from Oj 
Phys. Rev. Lett, 61:2546-2549, Nov 1988. 

[5] Serge Haroche and Daniel Kleppner. Gavity Quantum Electrodynamics. Physics To¬ 
day, 42(l):24-30, January 1989. 

[6] Eli Yablonovitch. Inhibited Spontaneous Emission in Solid-State Physics and Electronics, 

Phys. Rev. Lett., 58:2059-2062, May 1987. 


12 














[7] Peter Lodahl, A. Floris van Driel, Ivan S. Nikolaev, Arie Ir- 

man, Karin Overgaag, Daniel Vanmaekelbergh, and Willem L. Vos. 

Controlling the dynamics of spontaneous emission from quantum dots by photonic crystals. 
Nature, 430(7000):654-657,2004. 

[8] E A. Inam, T. Gaebel, C. Bradac, L. Stewart, M. J. 

Withford, J. M. Dawes, J. R. Rabeau, and M. J. Steel. 

Modification of spontaneous emission from nanodiamond colour centres on a structured surface, 
New Journal of Physics, 13(7):073012,2011. 

[9] Dmitri Toptygin. Effects of the Solvent Refractive Index and Its Dispersion on the Radiative Deca} 
Journal of Fluorescence, 13(3):201-219, 2003. 

[10] G. Nienhuis and C.Th.J. Alkemade. Atomic radiative transition probabilities in a continuous medi 
Physica B+C, 81(1):181 - 188,1976. 

[11] Roy J. Glauber and M. Lewenstein. Quantum optics of dielectric media, Phys. Rev. 

A, 43:467-491, Jan 1991. 

[12] Stephen M. Barnett, Bruno Huttner, and Rodney Loudon. 

Spontaneous emission in absorbing dielectric media, Phys. Rev. Lett., 68:3698- 
3701, Jun 1992. 

[13] Stephen M Barnett, Bruno Huttner, Rodney Loudon, and Reza Matloob. 

Decay of excited atoms in absorbing dielectrics. Journal of Physics B: Atomic, Molec¬ 
ular and Optical Physics, 29(16):3763,1996. 

[14] A. M. Stoneham. Theory of Defects in Solids. Glarendon Press:Oxford, 1985. 

[15] Mark Fox. Optical Properties of Solids, 2nd Ed. OUP:Oxford, 2012. 

[16] Pedro de Vries and Ad Lagendijk. Resonant Scattering and Spontaneous Emission in Dielectrics: R 
Phys. Rev. Lett., 81:1381-1384, Aug 1998. 


13 



[17] Michael Fleischhauer. Spontaneous emission and level shifts in absorbing disordered dielectrics ai 
Phys. Rev. A, 60:2534-2539, Sep 1999. 

[18] Michael E. Crenshaw and Charles M. Bowden. 

Effects of Local Fields on Spontaneous Emission in Dielectric Media, Phys. Rev. 

Lett, 85:1851-1854, Aug 2000. 

[19] P. R. Berman and P. W. Miloimi. Microscopic Theory of Modified Spontaneous Emission in a Dielec 
Phys. Rev. Lett., 92:053601, Feb 2004. 

[20] Hao Fu and P. R. Berman. Microscopic theory of spontaneous decay in a dielectric, 

Phys. Rev. A, 72:022104, Aug 2005. 

[21] D.V. Kuznetsov, Vl.K. Roerich, and M.G. Gladush. 

Local field and radiative relaxation rate in a dielectric medium, Journal of Ex¬ 
perimental and Theoretical Physics, 113(4):647-658,2011. 

[22] G. L. J. A. Rikken and Y. A. R. R. Kessener. 

Local Field Effects and Electric and Magnetic Dipole Transitions in Dielectrics. 

Phys. Rev. Lett., 74:880-883, Feb 1995. 

[23] P. Lavallard, M. Rosenbauer, and T. Gacoin. Influence of surrounding dielectrics on the spontaneoi 
Phys. Rev. A, 54:5450-5453, Dec 1996. 

[24] G. Lamouche, P. Lavallard, and T. Gacoin. Optical properties of dye molecules as a function of the 
Phys. Rev. A, 59:4668-4674, Jun 1999. 

[25] Frank J. P. Schuurmans, Dermis T. N. de Lang, Ger¬ 
ard H. Wegdam, Rudolf Sprik, and Ad Lagendijk. 

Local-Field Effects on Spontaneous Emission in a Dense Supercritical Gas. Phys. 

Rev. Lett., 80:5077-5080, Jun 1998. 


14 



[26] R. S. Meltzer, S. R Feofilov, B. Tissue, and H. B. Yuan. 

Dependence of fluorescence lifetimes of Y 2 O 3 : Eu^'*' nanoparticles on the surrounding medium. 
Phys. Rev. B, 60:R14012-R14015, Nov 1999. 

[27] G. Manoj Kumar, D. Narayana Rao, and G. S. Agarwal. 

Measurement of Local Field Effects of the Host on the Lifetimes of Embedded Emitters. 

Phys. Rev. Lett., 91:203903, Nov 2003. 

[28] G. Manoj Kumar and D. Narayana Rao. Modification of the spontaneous emission lifetime of Tb^"'' 
Optical Materials, 31 (9): 1343-1345,2009. 

[29] Ghang-Kui Duan, Hongli Wen, and Peter A. Farmer. 

Local-field effect on the spontaneous radiative emission rate, Phys. Rev. B, 

83:245123, Jun 2011. 

[30] Sander F. Wuister, Gelso de Mello Doneg, and Andries Meijertnk. 

Local-field effects on the spontaneous emission rate of GdTe and GdSe quantum dots in dielectric 
The Journal of Chemical Physics, 121(9):4310-4315,2004. 

[31] Ghang-Kui Duan and Michael F. Reid. Local field effects on the radiative lifetimes of Ge^’*' in differ 
Current Applied Physics, 6(3):348-350,2006. 

[32] Aleksander Zych, Matthijs de Lange, Gelso de Mello Doneg, and Andries Mei- 
jerink. Analysis of the radiative lifetime of Pr^’*’ d-f emission, Journal of Applied 
Physics, 112(1):-, 2012. 

[33] Ksenia Dolgaleva and Robert W. Boyd. Local-field effects in nanostructured photonic materials. 
Adv. Opt. Photon., 4(l):l-77, Mar 2012. 

[34] Seng-Tiong Ho and Prem Kumar. Quantum optics in a dielectric: macroscopic electromagnetic-fiel 
/. Opt. Soc. Am. B, 10(9):1620-1636, Sep 1993. 


15 



[35] J. H. Van Vleck. THE INFLUENCE OF DIPOLE-DIPOLE COUPLING ON THE DIELECTRIC COE 
Ann. N. Y. Acad. Sci., 40:293,1940. 

[36] R. Loudon. The Quantum Theory of Light. Oxford University Press:Oxford, 2010. 

[37] J. D. Jackson. Classical Electrodynamics. Wiley, second edition, 1975. 


16 


